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Abstract 


SARS-coronavirus (CoV). is a zoonotic agent derived from rhinolophid bats, in which a pleth- 
ora of SARS-related, conspecific viral lineages exist. Whereas the variability of virulence 
among reservoir-borne viruses is unknown, it is generally assumed that the emergence of 
epidemic viruses from animal reservoirs requires human adaptation. To understand the 
influence of a viral factor in relation to interspecies spillover, we studied the papain-like pro- 
tease (PLP) of SARS-CoV. This key enzyme drives the early stages of infection as it cleaves 
the viral polyprotein, deubiquitinates viral and cellular proteins, and antagonizes the inter- 
feron (IFN) response. We identified a bat SARS-CoV PLP, which shared 86% amino acid 
identity with SARS-CoV PLP; and used reverse genetics to insert it into the SARS-CoV 
genome. The resulting virus replicated like SARS-CoV in Vero cells but was suppressed in 
IFN competent MA-104 (3.7-fold), Calu-3 (2.6-fold) and human airway epithelial cells (10.3- 
fold). Using ectopically-expressed PLP variants as well as full SARS-CoV infectious clones 
chimerized for PLP, we found that a protease-independent, anti-IFN function exists in 
SARS-CoV, but not ina SARS-related, bat-borne virus. This PLP-mediated anti-IFN differ- 
ence was seen in primate, human as well as bat cells, thus independent of the host context. 
The results of this study revealed that coronavirus PLP confers a variable virulence trait 
among members of the species SARS-CoV, and that a SARS-CoV lineage with virulent 
PLPs may have pre-existed in the reservoir before onset of the epidemic. 
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Author summary 


Novel detection and sequencing technologies have greatly improved our knowledge of 
virus diversity in nature. Metaviromic screening of zoonotic animal reservoirs has become 
an established approach in pathogen surveillance and pandemic preparedness research. 
However, knowledge of viral genomes and host-associated viromes is insufficient to pre- 
dict zoonotic spillover events of reservoir-borne viruses. Phenotypic characterization of 
important viral functions will be necessary to identify virulence traits that determine the 
potential of viral emergence. As proof-of-principle, the present study demonstrates rele- 
vant functional differences between one of the main host immune antagonists of bat- 
borne viruses that belong to the same virus species as the epidemic agent of SARS. The 
antagonist, the papain-like protease, shows double action against IFN-mediated antiviral 
effects in epidemic SARS-coronavirus (binding and processing of ubiquitin), whereas the 
homologous protein in bat-borne viruses has only the processing function. This finding is 
surprising as the papain-like protease is a highly conserved protein domain that was not 
expected to vary among conspecific coronaviruses. PLP function may represent a variable 
virulence trait among reservoir-borne viruses. The preformed virulence of the primordial 
genetic lineage may have supported the emergence of SARS-CoV as a human epidemic 
agent. 


Introduction 


The coronaviruses (CoV; family Coronaviridae) are among the most relevant groups of viruses 
with zoonotic potential. CoVs are large, positive-sense, single-stranded RNA viruses that cause 
acute and prolonged infections ina variety of mammals and birds. Pathogenic human CoVs 
include members of the genus Alphacoronavirus, termed human coronavirus (HCoV)-NL63 
and HCoV-229E, as well as members of the genus Betacoronavirus, termed HCoV-OC43 and 
HCoV-HKUL. These endemic viruses cause upper and lower respiratory tract infections in 
humans worldwide. Past zoonotic descent can be inferred for HCoV-OC43 and -229E, respec- 
tively [1-4]. Actual zoonotic acquisition is known for two betacoronaviruses that both cause 
severe lung disease in humans. The Middle East respiratory syndrome coronavirus (MERS-- 
CoV) is a zoonotic agent that is frequently and repeatedly acquired by humans upon contact 
with dromedary camels in the Arabian Peninsula and parts of Africa [5, 6]. This virus seems to 
cause only limited human-to-human transmission, but is considered a major threat to global 
public health due to recurring nosocomial outbreaks that may facilitate onward adaptation to 
humans [7-9]. The severe acute respiratory syndrome (SARS)-CoV caused an epidemic with 
sustained human-to-human transmission during 2002 to 2003 in China and other countries, 
involving more than 8,000 notified infections with a case fatality proportion of about 10% [10, 
11]. 

Human SARS cases were derived from at least two independent zoonotic transmissions 
from a putative intermediary reservoir in feral carnivores. The majority of cases were part of 
one continuous chain of human-to-human transmission [12, 13]. Bats are now known to har- 
bor SARS-CoV strains that can directly infect primate cells [14-18]. SARS-CoV has become a 
paradigmatic subject to study pre-pandemic processes, both on an ecological and molecular 
level [14, 19-21]. 

Due to the zoonotic nature of many CoVs, there is an increasing interest to understand 
functional differences in the relative ways in which viruses deal with host cell defenses across 
different host species, particularly the innate immune response mediated by type I interferons 
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[22]. The type I interferon (IFN) response is an effective antiviral barrier that may limit zoo- 
notic cross-host infection in general terms [23, 24]. CoV infection is sensed by melanoma dif- 
ferentiation antigen 5 (MDAS5) and signaled via mitochondrial antiviral-signaling protein 
(MAVS), stimulator of IFN genes (STING), and IFN regulatory factor 3 (IRF-3), eventually 
leading to type I IFN gene transcription [25, 26]. 

Several IFN antagonist functions of zoonotic SARS-CoV are known. Within the viral struc- 
tural and accessory proteins, proteins 3b, 6, as well as the nucleocapsid proteim have been dem- 
onstrated to antagonize type I IFN (e.g., [27]). However, these viral proteins are expressed only 
after polyprotein processing, nonstructural gene expression, and subgenomic RNA transcrip- 
tion and, therefore, antagonize the downstream effects of IFN receptor signaling rather than 
the induction of IFN. Distinct from the above SARS-CoV IFN antagonists, the papain-like 
protease (PLP) is an IFN antagonist that constitutes a domain of the replicase polyprotein and, 
therefore, may be active at an early stage of the replication cycle to antagonize an upstream 
step of IFN induction. Additionally, unlike accessory proteins, which can vary greatly between 
CoV species, maintenance of PLP catalytic activity is critical to viral replication and is there- 
fore conserved across all CoVs [28, 29]. 

The coronavirus PLP proteins are multifunctional and encode a catalytic triad domain that 
catalyzes site-specific peptide cleavage of the viral polyprotein and the removal of both ubiqui- 
tin and IFN stimulated gene (ISG) 15.post-translational modifications. PLP protease activity 
catalyzes the processing of the replicase polyprotein at cleavage sites between nsp1/nsp2, nsp2/ 
nsp3, and nsp3/nsp4. PLP deubiquitinating (DUB) activity has been demonstrated in several 
CoV species, and acts directly andindirectly on several signal molecules in the IRF-3-depen- 
dent IFN induction pathway including retinoic acid inducible gene-I (RIG-I), tumor necrosis 
factor receptor-associated factor 3 (TRAF3), TANK-binding kinase 1 (TBK1) and STING [30, 
31]. K63-linked:polyubiquitin chains play a general role in signaling cascades of the proinflam- 
matory and IFN systems. K48-linked polyubiquitins label proteins for degradation by the pro- 
teasome and activate proinflammatory and antiviral factors. For instance, NF«B is activated by 
proteasomal degradation of its inhibiting factor IxB [32]. PLP DUB function also involves 
deISGylating activity, causing the removal of ISG15 modifications from viral and host proteins 
[33,34]. ISG15 is an IFN-inducible, antiviral protein that structurally resembles K48-linked 
di-ubiquitins. It can also be deconjugated by PLPs of other RNA viruses, in particular the 
ovary tumor (OTU) domain in the papain-like protease 2 (PLP2) of arteriviruses, and the L- 
gene-encoded OTU domain of nairoviruses [35, 36]. The protease recognition sequence 
LXGG is common to cleavage sites in the viral protein as well as ubiquitin and its derivatives. 
The DUB- and deISGylating activities inCoV PLPs should therefore be widely conserved. 

Due to the importance of ubiquitin-based innate immune functions, PLP functions may 
constitute a relevant predictor of the capability of reservoir-borne CoVs to overcome species 
barriers. PLP activity profiles may differ between relevant zoonotic CoV species. For instance, 
SARS-CoV has better ability to deconjugate K48-, as opposed to K63-linked polyubiquitins, 
whereas these activities are balanced in MERS-CoV [37-39]. The processing of ISG15 and 
K48-linked di-ubiquitin is more effective for SARS- than MERS-CoV PLP [37]. The blocking 
of induction of IFN by DUB activity was also confirmed for MERS-CoV, but in contrast to 
SARS-CoV, this inhibition is not independent of PLP’s protease activity [40, 41]. Interestingly, 
MERS-CoV is more sensitive to the effects of IFN than SARS-CoV [42, 43]. 

Unfortunately, the PLP activity profile cannot be derived from phylogenetic relatedness. 
For instance, above-mentioned studies found the distantly related HCoV-NL63 and SARS-- 
CoV to be similar in essential features such as protease-independent, DUB-mediated IFN 
antagonism, while the MERS-CoV that is much closer related to SARS-CoV only inhibits IFN 
induction when the protease function is intact (e.g., [34, 37]). Direct studies of PLP functions 
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of reservoir-borne viruses are therefore necessary to help us determine if there are differences 
in intrinsic virulence or virus-host interactions. 

In view of the complexity of PLPs interactions with innate immunity, functional studies 
have to take the whole viral replication cycle into account. To date, the DUB functions of 
SARS-CoV PLP have not been studied in the context of a replicating virus. Moreover, no stud- 
ies have so far compared PLP functions between members of one same viral species including 
natural variants existing in the zoonotic reservoir. Differences between reservoir-borne and 
epidemic viruses may uncover mechanisms that aid viral emergence of potentially pandemic 
strains. 

Based on epidemic and reservoir-borne variants of the species SARS-CoV, here we exem- 
plify functional differences in PLP domains. By reverse genetics, we show that the PLP of the 
epidemic SARS-CoV has an enhanced IFN antagonist function that is independent of PLP 
protease activity, and that is not present in the PLP of a bat-associated SARS-CoV. Additional 
mutagenesis studies in replicating virus context associate the functional difference to a more 
efficient binding of ubiquitin or ubiquitin-like modifiers. Against the general assumption that 
reservoir-associated viruses are highly adapted to their hosts, we find the PLP of the human 
epidemic virus to counteract IFN better than the bat-derived PLP even in bat cells. PLP func- 
tion is a viral virulence trait that varies among reservoir-borne viruses. 


Results 


In an earlier study we have described SARS-related CoVs in European (Bulgarian) bat species 
belonging to the genus Rhinolophus [44]. Fig 1A shows.a phylogeny of SARS-related beta- 
CoVs based on the PLP gene (981 bp fragment, genome position 4885 to 5829 in GenBank 
accession number AY310120). Based on standing classification criteria, the European bat- 
derived CoVs are conspecific with human SARS-CoV and in sister relationship to all Asian 
SARS-related CoVs. 

In addition, closely related viruses that were not conspecific with SARS-CoV but represent 
the closest phylogenetic outgroup to the species SARS-related CoV were discovered in Ghana- 
ian Hipposideros bats [45]. Hipposideros represents a sister genus to the typical SARS-CoV host 
Rhinolophus (Fig 1A). The PLP of human SARS-CoV is henceforth referred to as SA-PLP; the 
PLP of the conspecific European bat virus as SR-PLP (for SARS-Related); and the PLP of the 
sister species virus as SO-PLP (for SARS Outgroup). 


Comparison of PLP sequences 


An amino acid sequence alignment of the PLP region shows obvious similarities between 
SA-PLP and SR-PLP, and less so between these PLPs and SO-PLP. The PLP core domains in 
SA-PLP and SR-PLP each comprise 315 amino acids, and in SO-PLP 320 amino acids. SA-PLP 
and SR-PLP are 86% (271/315 aminovacids) identical. SO-PLP share 39% (125/324 positions 
including insertions/deletions) and 36% (118/324 positions including insertions/deletions) 
identical amino acids with SA- and SR-PLP, respectively (Table 1). A catalytic triad consisting 
of the three residues cysteine C1651, histidine H1812 and aspartic acid D1826 was previously 
shown to be responsible for cleavage of the SARS-CoV replicase polyprotein, and is present in 
SR-PLP (Fig 1B, grey arrows) [46]. In SO-PLP the aspartic acids (D1826) are replaced by an 
asparagine (N). This alternative type of catalytic domain was previously described for other 
cysteine proteases [47]. Another indispensable feature of SA-PLP is the zinc-binding domain, 
comprised of four cysteine residues, which connect the left- and right-hand domains of the 
papain-like fold by a zinc atom [48]. All these residues are also present within SR- and SO-PLP 
amino acid sequences (Fig 1B, marked with asterisks). 
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Fig 1. Phylogenetic and sequence-based analysis of the SARS-related bat coronavirus papain-like protease 
(SR-PLP). (a) Phylogeny of SARS-related beta-CoVs in the PLP gene (981 bp fragment) within the nonstructural 
protein 3. PLP genes characterized in the study are colored in red. The right-hand column shows the species 
classification of the included virus clades according to the International Committee on Taxonomy of Viruses (ICTV). 
Phylogenetic trees of SARS-related betacoronaviruses (CoVs) were calculated by the Neighbor Joining algorithm in 
Geneious under the assumption of a Tamura-Nei genetic distance model. Symbols correspond to the respective host 
species (human, civet and bat). The scale bar refers to the genetic distance. The SARS-outlier CoV (SO-CoV) was 
identified in a Ghanaian Hipposideros bat. SO-CoV belongs to a novel unclassified beta-CoV species. HCoV: human 
CoV, FRA: SARS Frankfurt strain, BtCoV: bat CoV. The accession numbers are as follows: HCoV_SARS/FRA: 
AY310120, Civet CoV_SARS: AY572034, BtCoV_Rp3: DQ071615, BtCoV_Rm1: DQ022305, BtCoV_Bulgarian: 
GU190215, BtCoV_Ganaian: MG916963, HCoV_MERS/EMC: JX869059. (b) Amino acid sequence alignment for the 
comparison of SR-PLP to SA-PLP. The alignment is based on the amino acid codes by the Blosum62 algorithm in the 
Geneious 6 software package. The SO-CoV derived PLP (SO-PLP) was included as an outlier PLP. Yellow boxes 
indicate conserved residues in all sequences. The boxes in light grey indicate conserved residues in only two sequences. 
Residues that form the catalytic center are indicated by grey arrows below the sequences. The catalytic cysteine, which 
was mutated to alanine in the course of this study, is highlighted in red. The ubiquitin-binding methionine at amino 
acid position 209, which was mutated to arginine (M209R) in this study, is marked in blue. Zinc-binding residues, 
important for the three dimensional PLP structure, are indicated by asterisks above the sequences. C1651 numeration 
refers to the position in the SARS-CoV ppla already used before [46]. Residues framed in black indicate the binding 
sites of the inhibitor compound 3e, which was used in the course of this study. SA: SARS; SR: Bulgarian; SO-PLP: 
Ghanaian. 


https://doi.org/10.1371/journal.ppat.1007296.g001 


Comparison of protease activities 


To functionally compare the PLPs, protease activities were assessed by a trans-cleavage assay 
[49]. The assay was based on coexpression of a SARS-CoV nsp2/3-GFP substrate with the 
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Table 1. Amino acid identity and similarity matrix. 





PLP Identity or similarity % identity or similarity with: 








SA-PLP SR-PLP SO-PLP 





SA-PLP Identity 





100 


























Similarity 100 

SR-PLP Identity 86.03 100 
Similarity 94.60 100 

SO-PLP Identity 38.58 36.42 100 
Similarity 58.64 58.95 100 








The identities and similarities of the listed proteins is based on an amino acid alignment using the BLOSUM62 substitution matrix and a threshold of 1. 


https://doi.org/10.1371/journal.ppat.1007296.t001 


respective PLPs, testing the cleavage of substrate into truncated products nsp2 and nsp3-GFP 
[49]. During establishment of the assay we noticed a considerable level of mRNA splicing 
while expressing SR-PLP under the control of a chicken B-actin promoter (S1 Fig). Codon- 
optimized constructs were therefore generated. The truncated nsp3-GFP product was detect- 
able by Western blot using anti-GFP epitope tag antibodies with all PLPs (Fig 2A; lanes 2, 4 
and 6). To confirm that the protease activity of SR-PLP and SO-PLP depended on the same 
typical catalytic domain as in SA-PLP, the catalytic cysteines of each PLP (C1651) were 
changed to alanines. For SA-PLP this mutation was previously shown to abolish PLP activity 
[46]. The mutants are henceforth referred to as CA-mutants. As expected, each PLP CA- 
mutant was unableto process the SARS nsp2/3-GFP substrate (Fig 2A; lanes 3, 5 and 7). To 
confirm PLP expression in all cases, Western blots were done on the same cell lysates using 
antibodies directed against the FLAG epitope tag fused to PLPs. It was found that the expres- 
sion levels of all PLPs were equal (Fig 2A, lower panel). 

To enable a quantitative comparison of protease activities, a PLP biosensor luciferase assay 
was done in which a split Firefly luciferase is coexpressed together with the PLP of interest. 
Upon cleavage of an LXGG protease cleavage site in the split luciferase construct, luciferase 
activity is reconstituted and measured after equilibration with cell membrane-penetrating 
luciferase substrate [50]. A time-course experiment confirmed that all PLPs had similar prote- 
ase activities ranging between 1- to 6-fold within 2 to 6 h when compared to CA-mutants cor- 
responding to each PLP (Fig 2B). 

To investigate if the amino acid differences between the PLPs affect the stereostructure of 
the catalytic site, a small molecule competitive inhibitor known to be specific for the PLP cata- 
lytic site, named 3e [51], was tested side-by-side on the PLPs. Inhibition of protease activity by 
the inhibitor was successful for both SA-PLP and SR-PLP, indicating structural similarity of 
both catalytic sites. The inhibitor was slightly more efficient for SA-PLP (ECs9 = 26.20 uM) 
than for SR-PLP (ECs9 = 30.28 uM). which is plausible because the inhibitor was designed to 
target a beta-loop structure (BL2) of SA-PLP located close to the catalytic site of the protease 
(Fig 1B). The inhibitor had very low efficiency towards SO-PLP, whose inhibitor-binding site 
has only 46% amino acid identity (6/13 amino acids identical to SA-PLP) to the inhibitor target 
site in SA-PLP (Fig 2C). 


Comparison of DUB functions 


To compare DUB activities, HEK-293T cells were cotransfected with increasing doses of plas- 
mids encoding each PLP along with constant doses of plasmids encoding HA-tagged ubiqui- 
tin. The decrease of protein ubiquitination conferred by PLP was determined by Western blot 
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Fig 2. SR-PLP had conserved-protease cleavage-, DUB- and deISGylating activities. (a) HEK-293T cells were transfected with empty vector 
plasmid (EV) or plasmids expressing either wild type (WT)—or catalytic mutant (CA) PLPs. SARS-CoV nsp2/3-GFP was cotransfected 
simultaneously. Lysates were harvested at 16 hours post transfection (hpt), and gene expression was analyzed by Western blotting. (b) A 
biosensor assay was applied for the detection of PLP cleavage activity. Cells were cotransfected with pGlo Firefly luciferase, and either WT-PLP, 
CA or EV plasmids in 96-wells. At 14 hpi, cells were incubated with GloSensor reagent and luminescence was detected. (c) To investigate the 
activity spectrum of a SA-PLP protease inhibitor, one hour after the GloSensor incubation, 12.5, 25 or 50 uM of compound 3e or DMSO were 
added. PLP activities were analyzed in relation to the different amounts of compound 3e at 4 h post treatment. Values were normalized to the 
respective DMSO-treated WT-PLP. Biosensor assays were performed in triplicate and repeated three times independently. Error bars indicate 
standard deviations of the means. Statistical significance between DMSO and inhibitor-treated cells or cells transfected with the CA-PLPs, 
respectively, was determined using one-way ANOVA and Sidak post hoc test. Statistically significant differences are indicated by asterisks (p> 
0.05 not significant (ns), p< 0.05 significant (*), p< 0.01 very significant (**), p< 0.001 highly significant (***)). (d) HEK-293T cells were 
transfected with EV or plasmids expressing either WT- or CA-PLPs. For analysis of DUB activity WT-PLP plasmids were transfected in 
increasing amounts of plasmids (50 ng, 100 ng and 200 ng per 12-well). Control plasmids (200 ng) EV and CA-PLP were transfected for 
comparison. HA-ubiquitin (HA-ub) was coexpressed in all samples. Lysates were harvested at 18 hpt, and gene expression was analyzed by 
Western blotting. (e) For analysis of deISGylating activity pISG15-myc and the conjugation machinery (UbcH8, Ube1L, and Herc5) of ISG15 
were coexpressed in all samples. Lysates were harvested at 18 hpt, and gene expression was analyzed by Western blotting. Western blot 
experiments were repeated three times independently and one representative blot is shown. B-Actin was applied as loading control. 


https://doitorg/10.1371/journal.ppat.1007296.g002 





using anti-HA antibodies. Each PLP deconjugated ubiquitin in a dose-dependent manner sug- 
gesting that the PLPs have comparable DUB efficiencies (Fig 2D). 

Because it is known that SA-PLP has deISGylating activity [34] the efficiency to deconjugate 
ISG15 from cellular proteins was determined by deISGylation assay. The PLPs were coex- 
pressed with myc-tagged ISG15 in HEK-293T cells. The extent of deISGylated proteins was 
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determined by Western blot using anti-myc antibodies. Each PLP was highly efficient in 
deconjugating ISG15 from the cellular proteins (Fig 2E). 

Taken together, these results suggested comparable levels of protease activity and identical 
DUB and deISGylating activities of SA- and SR-PLP when overexpressed in ahuman cell con- 
text. Differential efficiency towards a PLP inhibitor hint at structural differences even between 
PLPs from conspecific viruses that occur in zoonotic reservoirs. However, complete failure of 
the inhibitor was only seen with SO-PLP that falls outside the limits of current CoV species 
classification. 


Protease function but not IFN antagonism is equivalent for SARS-CoV 
PLPs in full virus context 


To quantitatively compare the functions of the two closely related PLPs (SA- and SR-PLP) in 
the context of the full virus replication cycle, we constructed a chimeric SARS-CoV in which 
SA-PLP was replaced by SR-PLP (Fig 3A). Viral plaques with similar morphologies were 
observed for the recombinant wild type virus (rSCV, Fig 3B) and the chimeric virus (SR-PLP- 
rSCV, Fig 3C), indicating that the SR-PLP wasable to functionally compensate the SA-PLP. 
As already suggested by protease cleavage assays, SR-PLP-rSCV was slightly less sensitive 
against protease inhibitor 3e than SA-PLP-rSCV (rSCV: ICs9 = 2.36 uM; SR-PLP-rSCV: ICso = 
11.02 uM), providing additional evidence for functional and structural integrity of SR-PLP in 
the context of SARS-CoV replication (Fig 3D). To exclude the possibility that 3e was cytotoxic, 
a cell viability assay was conducted) The number of viable cells decreased with increasing 
amounts of inhibitor to a minimum of 80% at the highest dose of 50 uM (S2 Fig). This con- 
firms the specific action of 3e. 

To obtain amore quantitative comparison of SA- and SR-PLP during virus replication, 
multistep growth curve experiments were done for both viruses in Vero cells. Both viruses 
grew to the same titers at.all tested time points (range: 8.2x10° PFU/ml to 8.2x10° PFU/ml; 8, 
14, 24, 48 hours post infection [hpi], Fig 4A). Notably, growth curves differed when both 
viruses were grown in the type I IFN-competent primate cell line MA-104. SR-PLP-rSCV grew 
to significantly lower titers than rSCV (general linear regression model, p = 0.038/R- 
square = 0,612; differences 3.7-fold at 14 hpi and 3.6-fold at 24 hpi, Fig 4B). The reduced 
growth of SR-PLP-rSCV in MA-104 cells may indicate a less efficient viral counteraction 
against type I IFN. 

In order to compare the IFN sensitivity of both viruses, Vero cells were treated with a 
defined concentration of pan-species IFN-c. A significantly increased IFN sensitivity (4.2-fold 
at 100 IU/ml; p = 0.008 in a two-sided t test) compared to rSCV confirmed that the anti-IFN 
activity of SR-PLP was decreased in primate cells (Fig 4C). 

To further confirm that growth of SR-PLP-rSCV is also reduced in context of the human 
cell environment, we performed multistep growth curve experiments in the type I IFN compe- 
tent human lung epithelial cell line Calu-3 (Fig 4D). Again, SR-PLP-rSCV grew to significantly 
lower titers than rSCV (2.6 fold at 24 hpi [p = 0.023] and 1.5 fold at 48 hpi [p = 0.016] in a two- 
sided t test). A 1.5-fold increased detection of IFN-8 mRNA expression levels in SR-PLP- 
rSCV- compared to rSCV-infected Calu-3 cells further confirmed that the SR-PLP is less effi- 
cient in blocking IFN induction (Fig 4E). 

To better reflect the human respiratory tract and to generalize the notion that SR-PLP- 
rSCV grows less efficiently in presence of an active type I IFN response, multistep growth 
curve experiments were conducted in human airway epithelial cells (HAE). In accordance 
with our previous studies [42] virus growth was generally delayed in HAE compared to the 
monoclonal primate and human cell cultures. Importantly, a significantly reduced growth was 
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Fig 3. SR-PLP supported viral replication in the context of a recombinant; chimeric SARS-CoV (rSCV). (a) A chimeric rSCV containing 
SR-PLP (SR-PLP-rSCV) was generated by reverse genetics to investigate SR-PLP functions in the context of a replicating virus. SR-PLP 
(purple) was inserted at the genomic position of SA-PLP (blue). The plaque morphology was analyzed in Vero E6 cells. Therefore, cells were 
infected with (b) rSCV and.(c), SR-PLP-rSCV (MOI 0.01) and overlaid with a highly viscous medium. At 3 dpi, cells were fixed and stained 
with crystal violet. (d) To investigate structural integrity of the SR-PLP domain in the molecular context of the SARS-CoV nonstructural 
protein 3, cells were infected with either rSCV or SR-PLP-rSCV (MOI 0.0001) and treated with the SA-PLP protease inhibitor 3e (24-well 
format). After 1h, cells were washed twice with PBS and either DMSO (01M) or in DMEM serially diluted compound 3e (3.125, 6.25, 12.5, 25 
and 50 uM) was added. Supernatants were collected at 24 hpi. For virus quantification a real-time RT-PCR for genomic SARS-CoV RNA was 
performed. The experiment was done in triplicate and repeated twice. Error bars indicate the standard deviations of the means. Statistical 
significance between DMSO- and compound-treated cells was determined using one-way ANOVA and Sidak- or Games-Howell post hoc 
tests for rSCV and SR-PLP-rSCV, respectively. The 50% effective concentrations (ECso) were rSCV: ICs = 2.36 uM and SR-PLP-rSCV: ICs9 = 
11.02 uM. 


https://doi.org/10.1371/journal.ppat.1007296.g003 





detected for SR-PLP-rSCV compared to rSCV at 96 hpi (10.3 fold; p = 0.006 in a two-sided t 
test) providing further confirmation that the anti-IFN activity of SR-PLP may be decreased in 
type I IFN competent cells (Fig 4F). 


Analysis of differential sensitivity to IFN 


The treatment with IFN-o in the above experiment (Fig 4C) would broadly affect IFN induc- 
tion, signaling, and response. Because PLP IFN antagonism functions have been linked to IRF- 
3 function, a sensitive assay for IRF-3 nuclear translocation was established. Nuclear transloca- 
tion of an overexpressed IRF-3/GFP fusion protein was stimulated by superinfection with Rift 
Valley fever virus clone 13 (RVFV Cl 13), an RVFV-mutant devoid of the IFN induction antag- 
onist NSs. RVFV CIl 13 is known to trigger a strong IFN response [52, 53]. The proportion of 
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Fig 4. Chimeric SR-PLP-rSCV grew less efficient in presence of type I interferon (IFN). Virus growth was compared in (a) type I IFN- 
deficient (Vero) primate cells, and (b) IFN-competent (MA-104) primate cells. For virus growth kinetics, cells were infected with rSCV and 
SR-PLP-rSCV (MOI 0.01). Supernatants were taken at 0, 8, 14, 24 and 48 hpi, and viral replication was determined by a plaque titration assay. 
Growth experiments in Vero and MA-104 cells were done in triplicate and repeated twice. Error bars indicate the standard deviations of the 
means. Infectious particle production of rSCV and SR-PLP-rSCV was compared-using SPSS Version 23.0.0.0 and a general linear regression 
model. Growth of both viruses did not significantly differ in Vero cells (p= 0.929 and R-square = 0.670). In MA-104 cells growth of the viruses 
significantly differed (p = 0.038 and R-square = 0.612). (c) Vero cells were treated with 100 IU/ml of recombinant pan-species IFN-a. At 16 h 
post IFN treatment, cells were infected with rSCV and SR-PLP-rSCV (MOI 0.01), respectively. Supernatants were taken at 24 hpi, and viral 
replication was determined by a plaque titration assay. The experiment was performed in triplicate and repeated twice. Error bars indicate the 
standard deviations of the means. One representative experiment is shown. (d) Virus growth was compared in IFN-competent human cells 
(Calu-3) as described above. The experiment was performed in triplicate. (e) To determine IFN-B expression, Calu-3 cells were infected with 
rSCV, SR-PLP-rSCV or IFN-inducing RVFV Cl 13 (control of IEN-B expression) at an MOI of 1. Total mRNA was extracted from cell lysates at 
24 hpi. IFN-B expression was determined using quantitative real-time PCR analysis] The mean fold change in IFN-B expression was calculated 
using TATA-bindi ng protein (TBP) expression as a reference gene and the 2-“““' analysis method [55]. The experiment was done in 
quadruplicates. Statistical significance between the indicated groups was determined using a two-sided t test. (f) Human airway epithelial cells 
(HAE) were infected with rsCV and SR-PLP-rSCV with an absolute infectious dose of 40,000 PFU. At 0, 48, 72 and 96 hpi samples were taken 
from and viral replication was determined by a plaque titration assay. The experiment was done in duplicate and repeated three times 
independently. Statistical significance in (d-f) was determined using a two-sided t test. 


https://doi.org/10.1371/journal.ppat.1007296.g004 


cells with nuclear translocation of GFP signal was counted microscopically. As summarized in 
Fig 5A and 5B, ectopic expression of SA- and SR-PLP blocked the nuclear translocation of 
IRF-3 to comparable levels (SA-PLP: 18% and SR-PLP: 20% of cellular IRF-3/GFP fusion pro- 
teins located in the nucleus). The SO-PLP of the outlying virus blocked the nuclear transloca- 
tion of IRF-3 even more efficiently (3% of cellular IRF-3/GFP fusion proteins located in the 
nucleus). For all PLPs, CA-mutants were included in the experiment to determine whether 
anti-IFN effects depended on protease function. The inhibitory capacity of the SA- and 
SO-PLP CA-mutants were strongly reduced, but still detectable at significant levels (SA-PLP: 
71% and SO-PLP: 45% of cells with IRF-3/GFP fusion protein located in the nucleus). In 
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Fig 5. Protease-dependent and protease-independent IFN-antagonistic functions of PLPs. (a) Primate cells (MA-104) were transfected 
with plasmids expressing 1 ug of WT- or CA-PLPs. EV was applied as a control. GFP-IRF-3 was coexpressed in each sample. At 24 hpi, cells 
were infected with IFN-inducing, recombinant Rift Valley Fever virus clone 13 (RVFV Cl 13) at an MOI of 5. Cells were fixed at 8 hpi. 
Immunofluorescence pictures show representative results. White arrows indicate representative cells. (b) For quantification of IRF3 
translocation events at least 4 representative pictures were taken. The total number of cells positive for GFP-IRF-3 and the respective PLP were 
determined and the number of GFP-IRF-3 nuclear translocation was calculated. The result represents one of three independently performed 
experiments. (c, d, e) An IFN-f promoter activation assay was conducted in HEK-293T cells to determine the anti-IFN activities of (c) SA-PLP 
and (d) SR-PLP and (e) SO-PLP. Cells were transfected with plasmids expressing MDAS (IFN stimulator), Firefly- (detection of IFN-B 
promoter activity) and Renilla- (detection of the general transcription level) luciferases. 50 ng/24-well of EV were transfected and samples 
were applied as induction controls. WT- and CA-PLPs were coexpressed in a dose-dependent manner (10, 25 and 50 ng/24-well). At 17 hpt, 
cells were lysed and the luciferase activity was measured. Results are presented as induction relative to EV. The graph shows results of three 
independently performed experiments, which were all done in triplicate. Expression of the PLPs was confirmed by Western blotting. B-Actin 
was applied as a loading control. The error bars in each graph indicate the standard deviations of the means. Statistical significance between 
infected/induced EV samples and PLP-expressing cells was determined using one-way ANOVA and Dunnett-T3 post hoc test. 


https://doi.org/10.1371/journal.ppat.1007296.g005 





contrast, the SR-PLP CA-mutant had essentially lost its ability to block IRF-3 nuclear translo- 
cation (94% of cells with nuclear IRF-3/GFP signal). 
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To see whether the PLP activity against IRF-3 translocation determined IFN counteraction, 
the abilities to block IFN promoter activation were compared in a reporter gene assay. All 
PLPs efficiently blocked IFN-B promoter activation in a dose-dependent manner (6.3%, 2.3% 
and 2.2% residual IFN-B promoter activation, respectively, at plasmid concentrations of 50 ng 
for SA-PLP, SR-PLP and SO-PLP; Fig 5C-5E). However, differences were observed with CA- 
mutants. The SA-PLP and SO-PLP CA-mutants still blocked the IFN-B promoter activity in a 
dose-dependent manner, albeit with reduced efficiency as compared to wild type (SA-PLP: 
26.4% [Fig 5C] and SO-PLP: 12.9% [Fig 5E] residual IFN-B promoter activation at a plasmid 
concentration of 50 ng). The SR-PLP CA-mutant lost any significant blocking functions 
against IFN-B promoter activation (Fig 5D). 

The differential capabilities of SA-PLP, SR-PLP and SO-PLP to counteract IFN induction 
may be attributable to differential capabilities to inactivate ubiquitin or ubiquitin-like modifi- 
ers such as ISG15. SR-PLP may have to rely on protease-dependent cleavage, while SA-PLP 
and SO-PLP may be able to provide additional IFN-antagonistic activity by binding of ubiqui- 
tin or ubiquitin-like modifiers. 


SA-PLP’s protease-independent IFN-antagonistic activity depends on the 
interaction between SA-PLP and ubiquitin 


To understand if within the SARSr-CoV species the SA-PLP as compared to SR-PLP provides 
additional anti-IFN activity depending on ubiquitin-binding but not cleavage, we decided to 
introduce a modification in the ubiquitin-binding surface of SA-PLP that should not affect the 
protease-processing activity. Amino acid residue M209 of SA-PLP interacts directly with a 
binding patch.comprised of 144, V70, L8, R167 and D168 of ubiquitin based on co-crystalliza- 
tion studies [38, 54]. In the MERS-CoV PLP, replacement of methionine at a corresponding 
position by the bulkier arginine residue was found to prevent the interaction between ectopi- 
cally expressed PLP and ubiquitin [41]. We thus decided to study the effect of a mutation that 
selectively destabilizes ubiquitin-binding in SARS-CoV via an M209R mutation. 

In order to specifically focus on the consequences of the mutation on IFN antagonism, we 
had to exclude that M209R disturbs the protease activity of the PLP. Expression plasmids were 
generated, and protease activities were compared by trans-cleavage assays, confirming that the 
mutation had no impact onthe cleavage activity (Fig 6A). M209R PLP processed the nsp2/3 
substrate with the same efficiency as SA-PLP. 

To confirm that ubiquitin-binding of SA-PLP is reduced by the M209R mutation, a DUB 
assay was conducted with wild type and M209R SA-PLP expression plasmids as previously 
described. DUB efficiency of M209R PLP (Fig 6B, lane 4) was reduced compared to the wild 
type PLP (Fig 6B, lane2) as indicated by the increased number of ubiquitinated protein bands 
in lane 4. 

The reduced ubiquitin-binding efficiency should result in a decreased IFN-antagonistic 
activity of M209R PLP. We thus compared IFN-B promoter activation in a reporter gene assay. 
Both, the M209R and the corresponding M209R CA mutant, had less anti-IFN activity com- 
pared to SA-PLP (Fig 6C). The combination of M209R and CA (Fig 6C, right black column) 
resulted in the complete depletion of significant anti-IFN activity. 

To explore the consequences of the disturbed ubiquitin-binding on virus growth, the 
M209R mutation was introduced into rSCV. As expected, the resulting virus, termed 
rSCV-PLP/M209R, grew to titers similar to rSCV in a multistep growth curve in Vero cells 
(Fig 6D). As observed with SR-PLP, the congruence with wild type changed when both viruses 
were grown in type I IFN competent MA-104 cells. In MA-104 cells rSCV-PLP/M209R grew 
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Fig 6. The protease-independent IFN-antagonistic activity of SA-PLP is related to the interaction between SA-PLP and ubiquitin. (a) A 
trans-cleavage assay was conducted in HEK-293T cells to determine the catalytic activity of SA-PLP and the SA-PLP mutant M209R, which 
had the ubiquitin-binding residue methionine (M) at amino acid position 209 mutated to arginine (R). (b) For analysis of DUB activity, HEK- 
293T cells were transfected with HA-ub and 200 ng of plasmids expressing either SA-PLP, the SA-PLP mutant M209R or the respective CA- 
mutants. Lysates were harvested at 18 hpt, and gene expression was analyzed by Western blotting. (c) IFN-B promoter activity was determined 
using an IFN-B promoter activation assay. 50 ng of SA-PLP encoding plasmids were coexpressed with the luciferase genes and IFN-inducing 
MDAS. Statistical significance between infected/induced EV samples and PLP-expressing cells was determined using one-way ANOVA and 
Dunnett-T3 post hoc test. Virus growth of rSCV and the modified rsCV-PLP/M209R, carrying the M209R mutation, was compared in type I 
IFN-deficient primate cells (Vero, d), and IFN-competent primate cells (MA-104, e). For virus growth kinetics, cells were infected with rSCV 
and rSCV-PLP/M209R at an MOI of 0.01. Each experiment was performed in triplicate. MA-104 growth experiment was repeated twice. 
Particle production of rsCV and rSCV-PLP/M209R was compared using SPSS Version 23.0.0.0 and a general linear regression model. Growth 
of both viruses did not significantly differ in Vero cells (p = 0.613 and R-square = 0.020). In MA-104 cells growth of the viruses significantly 
differed (p = 0.034 and R-square = 0.217). (f) In order to determine the sensitivity of rsCV and rSCV-PLP/M209R towards type I IFN, Vero 
cells were treated with 100 IU/ml of recombinant pan-species IFN-a. At 16 h post treatment, cells were infected with rSCV or rSsCV-PLP/ 
M209R (MOI 0.01) for 24 h. The experiment was performed in triplicate. Statistical significance between the indicated groups was determined 
using a two-sided t test. (g) Virus growth was compared in human-derived, IFN-competent cells (Calu-3). Cells were infected with rSCV or 
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rSCV-PLP/M209R at an MOI of 0.01. The experiment was done in triplicate. Statistical significance was determined using a two-sided t test. 
(h) To determine IFN-B expression, Calu-3 cells were infected with rSCV, rSCV-PLP/M209R or IFN-inducing RVFV Cl 13 (control of IFN-B 
expression) at an MOI of 1. Total mRNA was extracted from cell lysates at 24 hpi. IFN-B expression was determined using quantitative real- 
time PCR analysis. The mean fold change in IFN-B expression was calculated using TATA-b inding protein (TBP) expression as.a reference gene 
and the 24“ analysis method [55]. The experiment was done in quadruplicates. Statistical significance between the indicated groups was 
determined using a two-sided t test. 


https://doi.org/10.1371/journal.ppat.1007296.g006 
































to significantly lower titers than rSCV (general linear regression model, p = 0.034/R- 
square = 0.217; differences 1.9-fold at 14 hpi and 3.5-fold at 24 hpi, Fig 6E). 

To confirm that the destabilizing mutation in the ubiquitin-binding domain caused 
increased IFN sensitivity, Vero cells were pretreated with 100 IU of pan-species IFN-o (Fig 
6F). Upon treatment, the growth difference between both viruses was augmented, confirming 
a significantly increased IFN sensitivity of rSCV-PLP/M209R (4.3-fold at 0 [U/ml [p = 0.023] 
and 5.2-fold at 100 IU/ml [p = 0.007] in a two-sided t test). 

To further confirm the type I IFN-dependent growth differences of rSCV and rSCV-PLP/ 
M209R, multistep virus growth kinetics were done in the human lung epithelial cell line Calu- 
3 (Fig 6G). Compared to wild type virus-infected cells, the M209R mutation led to significantly 
lower titers in rSCV-PLP/M209R-infected Calu-3 cells (8.3 fold at 24 hpi [p = 0.007] and 2.8 
fold at 48 hpi [p = 0.003] in a two-sided t test). 

In order to directly detect the virus-dependent effect on IFN-B gene expression, Calu-3 cells 
were infected with both viruses with an MOI of 1 for 24 h and IEN-B mRNA upregulation was 
detected by quantitative real time RT-PCR assay (Fig 6H). To determine the fold-induction of 
the IFN- mRNA, the 2-“““*t method was applied with TATA» binding protein (TBP) as house- 
keeping gene [55]. TherSCV-PLP/M209R- infected cells had a more pronounced IFN-B 
mRNA expression level (3.2-fold difference between rSCV and rsCV-PLP/M209R at 24 hpi; 

p = 0.031 ina two-sided t test) compared to rSCV. The type I IFN-inducing RVFV Cl 13 was 
included as an IFN-B mRNA expression control showing that the cells had a fully functional 
IFN response. 

The phenotypic similarity between SR-PLP-rSCV and rSCV-PLP/M209R suggests that 
SA-PLP, but not SR-PLP, provides additional IFN antagonism via binding of ubiquitin or ubi- 
quitin-like modifiers. 


Differential IFN counteraction in cells of the SARS-CoV natural host 


As ubiquitin (and its derivatives) are conserved in mammalian species, we expected that ubi- 
quitin-dependent anti-IFN activities should be independent of the host species. As all mem- 
bers of the species SARS-CoV are carried by bats of the genus Rhinolophus, IFN-competent 
Rhinolophus alcyone lung epithelial cells (RhiLu) were made transgenic for the human ACE2 
receptor (RhiLu-hACEZ2) by lentiviral transduction [56], and were infected with rSCV and 
SR-PLP-rSCV (Fig 7A). In multistepgrowth curves, SR-PLP-rSCV grew to significantly lower 
titers compared to rSCV (13-fold lower at 24 hpi [p = 0.016] and 18-fold lower at 48 hpi 

[p = 0.005], respectively). Decreased viral growth of SR-PLP-rSCV in RhiLu-hACE2 bat cells 
may indicate more efficient IFN counteraction by rSCV compared to SR-PLP-rSCV, as in 
IFN-competent MA-104 cells. Remarkably, the growth differences were even more prominent 
in the bat cells compared to IFN competent primate cells (13-fold in RhiLu-hACE2 versus 
3.6-fold in MA-104 cells at 24 hpi). 

In order to investigate if the growth differences in bat cells were reflected by differences in 
IFN suppression, a quantitative real-time-PCR assay specific for the Rhinolophus IFN-B gene 
was established. To determine the specific IFN-8 mRNA expression levels, RhiLu-hACE2 cells 
were infected with rsCV, SR-PLP-rSCV and, as IFN induction control, RVFV Cl 13 (MOI of 
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Fig 7. SR-PLP-rSCV grew less efficiently in IFN-competent, bat-derived lung cells. (a) For virus growth kinetics Rhinolophus 
alcyone lung cells, which expressed the SARS-CoV receptor human angiotensin converting enzyme (ACE)-2, were infected with 
rSCV and SR-PLP-rSCV at an MOI of 0.1. Supernatants were taken at 0, 24, 48. and 72 hpi, and viral replication was determined by a 
plaque titration assay. The experiment was performed in triplicate. Error bars indicate the standard deviations of the means. (b) 
RhiLu-hACE2 cells were infected with rSCV, SR-PLP-rSCV or IFN-inducing RVFV Cl 13 (control of IFN-B mRNA expression) at 
an MOI of 1. Total mRNA was extracted from cell lysates 24 hpi. IEN-f expression was determined using quantitative real-time PCR 
analysis. The mean fold change in the IFN-f expression was calculated using f-actin expression as a housekeeping/reference gene 
and the 24“ analysis method [55]. The experiment was done in quadruplicates and repeated twice. One out of two experiments is 
shown. Statistical significance between the indicated groups was determined using a two-sided t test. 








https://doi.org/10.1371/journal.ppat.1007296:g007 


1). Total RNA was extracted from cell lysates at 24 hpi. To determine the fold-induction of the 
IFN-B mRNA expression levels, the 2-““““ method was applied with B-actin as housekeeping 
gene [55]. In comparison to rSCV, SR-PLP-rSCV infection led to significantly higher IFN-B 
expression levels (6.6 fold at 24 hpi [p = 0.001] in a two-sided t test, Fig 7B). 

In summary, these results suggest that the increased IFN sensitivity of SR-PLP is indepen- 
dent of the host cell species: 


Discussion 


SARS-CoV variants exist across Europe in widespread Rhinolophus bat species. It is important 
to understand whether these viruses constitute a risk for human infection [57]. Host tropism is 
mainly thought to be determinedby the spike protein, but studies have shown that CoV popu- 
lations in natural reservoirs can contain a plethora of spike variants, including variants that 
can directly mediate entry into human cells [16, 20]. Because of widespread recombination, 
spike proteins can be exchanged between viral genetic lineages in the reservoir. The spike pro- 
tein may not sufficiently represent the variability of virulence traits in the reservoir. 

Here we studied a conserved viral function related to host interaction by focusing on PLP, a 
multifunctional protein that drives essential steps in the infection cycle including cleavage of 
the viral polyprotein as well as deubiquitination and deISGylation. While these latter functions 
interfere with molecules of the IFN pathway and cytokine production, it was unknown to what 
extent they contribute to the replication phenotype of SARS-CoV. 

Our initial comparisons of PLP functions of epidemic and reservoir-associated virus were 
based on ectopically-expressed protein, confirming that essential properties such as cleavage of 
the viral polyprotein and sensitivity to a known inhibitor were almost identical between epi- 
demic- and reservoir-derived PLP. However, overexpression assays cannot reflect the whole 
complexity of the viral life cycle including the timing and compartmentalization of functions. 
To detect more the subtle differences that we expected to occur between natural virus variants, 
we generated an infectious clone carrying the reservoir-derived PLP in an otherwise 
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unmodified SARS-CoV backbone. Growth properties of the chimeric virus in IFN-deficient 
Vero cells confirmed that the reservoir-derived PLP effectively cleaved the viral polyprotein 
and supported recombinant virus replication kinetics comparable to wild type SARS-CoV. 
Only when grown in cells with a fully active IFN system, wild type SARS-CoV replicated better 
than the PLP-chimeric virus. The small but significant growth differences ranged between 2.6 
(Calu-3) and 10.3-fold (HAE) depending on the time point post infection and the applied cell 
cultures. We assume that other viral proteins with IFN antagonistic functions, like protein 6 
and 3b, might have compensated for the decreased immune-modulating functions of SR-PLP 
and PLP/M209R preventing a more pronounced growth difference [27]. Importantly, we 
showed that the effect was increased when artificial IFN was added to the cell cultures and was 
additionally accompanied by a more efficient blocking of IFN-8 mRNA upregulation, obtain- 
ing further confirmation for the specificity of the PLP function within the IFN response. 

We showed by studies on catalytically-inactivated PLP constructs that the anti-IFN function 
of the reservoir-derived SR-PLP depended on the protease function while SARS-CoV PLP 
exerted an additional, protease-independent activity in IFN evasion. The likely mechanism for 
this additional activity has been identified by in-vitro studies to involve binding of ubiquitin or 
ubiquitin-like modifiers [34, 54, 58]. Here we confirmed the functional link between DUB 
functions and IFN antagonism for the first time in a replicating CoV. Because the cytopatho- 
genic nature of SARS-CoV infection prevented sensitive cell-based assays for DUB- and delIS- 
Gylating functions, we relied.on knowledge of structure-function relationships from a study 
on equine arteritis virus (EAV), an arterivirus related to CoVs [35]. In EAV PLP2, a ubiquitin- 
binding surface distinct from the protease active site was identified and mutagenized without 
causing a loss of protease processing. A definite link to IFN antagonism was established by 
engineering of a destabilizing mutation into the ubiquitin-binding surface in an EAV infec- 
tious clone, causing increased induction of IFN during replication. By introduction of a 
homologous mutation in our SARS-CoV infectious clone, we found an even clearer loss of 
IFN antagonism than with EAV. The choice of the mutation was informed by another study 
that defined the primary ubiquitin-binding surface in SARS-CoV PLP and found position 
M209 to critically interact with a conformational binding patch (144, V70, L8) on ubiquitin 
but.to not affect protease activity [54]. Concordantly, we found no signs of catalytic inactiva- 
tion of the M209 mutant in Vero cells where it replicates like wild type virus, and a pattern of 
fully protease-dependent IFN antagonism like that of the reservoir-derived PLP. We conclude 
that ubiquitin binding constitutes.an anti-IFN function of SARS-CoV PLP that is independent 
of the conserved protease function, is phenotypically relevant for replication level and immune 
evasion, and is variable among viral variants. 

The IFN-related effects observed in the present study seem to be relevant as they resemble 
effects observed upon deletion of the 2’-O-methyltransferase activity, another essential repli- 
case function that modifies IFN-recognition [59]. In that study, a D130A mutation in the 
nsp16 2’-O-methyltransferase domain replicated like wild type in Vero cells but caused about 
10-fold lower replication in an IFN-competent human airway epithelial cell line. When engi- 
neered in a mouse-adapted backbone, the mutant was cleared faster from the lungs of infected 
mice, and caused significantly less weight loss. The present study did not use in vivo models to 
evaluate the consequences of the introduced mutations on replication. Mouse experiments 
may not reflect the role of PLP in host switching because ISG15 (expectably the most impor- 
tant ubiquitin-like modifier processed by SARS-CoV PLP [38]) is not conserved between 
humans, bats, and mice [60]. However, we can derive from the existence of SR-PLP in bats 
that the ubiquitin-related function of PLP is not essential for SARS-CoV replication in the nat- 
ural host, and functional variants reflect natural diversity. 
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While it has been suggested that SARS-CoV may have gained virulence for humans during 
human-to-human passage and adaptation [12, 61], the differences in immune interaction 
observed in the present study are not dependent on the host cells used. Even in bat cells from 
the natural host of SARS-CoV, epidemic SARS-CoV has additional anti-IFN activity. Differen- 
tial virulence via PLP may therefore have pre-existed within viral reservoir, rather than having 
evolved by positive selection upon adaptation within the human host. Therefore, we postulate 
that PLPs are a variable virulence trait among members of the species SARS-CoV that exist in 
natural reservoirs. 

Based on the PLP structure, there is huge potential for natural sequence variants to influ- 
ence ubiquitin-binding. For instance, PLP has a second ubiquitin-binding site enabling a bi- 
dentate binding mechanism that may explain the PLP preference for di-ubiquitins and ISG15 
[54]. Crystallization of PLP in complex with K48-linked di-ubiquitin identified the structure of 
the second binding site (Ub2), which seems to be more critical for K48-linked di-ubiquitin- 
binding than the primary binding surface [38]. The primary ubiquitin-binding surface (Ub1) 
is more critical for ISG15 binding [38,60]. The SARS-CoV strains found in reservoirs in 
China and Europe show several variants at both ubiquitin-binding sites that could be studied 
for effects on replication level or other virulence traits (Fig 8A). 

In silico modeling of the ubiquitin-binding sites of SA- and SR-PLP shows that sequence 
variations may even cause structural changes in- and adjacent to the binding sites (Fig 8B and 
8C). The fact that, apart from our observed functional differences, there may also be structural 
changes between two closely related PLPs of a single virus species (the species “SARS-related 
CoV”) is particularly intriguing as such differences were previously only seen among CoVs 
that belong to different virus genera. The PLP of HCoV-NL63, for example, has similar prote- 
ase-independent anti-IFN functions as the PLP of SARS-CoV whereas MERS-CoV PLP anti- 
IFN activity clearly depends on protease function [34, 37]. The existence of functional diversity 
in the reservoir makes sense in light of viral evolution and emergence. For instance, the capac- 
ity to modulate processing and substrate binding in separate domains of PLP opens the possi- 
bility for reservoir-associated:CoVs to adjust fitness or virulence levels based on changes in 
host population structure. Increased virulence may confer increased transmission rates, but 
may.come at the cost of host population decline. Based on our data it seems that the reservoir 
contains aconsiderable degree of functional variability that is based on highly conserved host 
functions such as the ubiquitin systems. Virus adaptation to the human host may not only be 
an evolved property at the beginning of virus emergence. Rather, the reservoir may contain 
pre-formed virulence traits that can be predicted by experiments informed by detailed knowl- 
edge of molecular virus functions, enabling new approaches to forecast emergence potential. 


Materials and methods 
Ethics statement 


Fecal pellets of Hipposideros bats were provided by author Christian Drosten. All fecal samples 
were collected non-invasively by trained field biologists and stored at -80°C until analysis. In 
detail, bats were caught with mist nets, which were checked at intervals of 5 min. Captured 
bats were freed from nets immediately and put into cotton bags for several minutes. While 
being kept in bags, bats produced fecal pellets that were collected and transferred into RNAla- 
ter RNA stabilization solution (QIAGEN, Hilden, Germany). Procedures were previously 
described in [62], and were consistent with guidelines of the American society of mammalo- 
gists [63] for the use of wild mammals in research and national guidelines for the capture, han- 
dling, and care of bats. For all capturing, sampling and sample export, permission was 
obtained from the Wildlife Division of the Ministry of Lands, Forestry, and Mines in Ghana 
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Fig 8. Variability in the PLPs of SARS-related CoVs found in the reservoirs in China and Europe. (a) Schematic representation of a PLP 
amino acid sequence alignment (315 amino acids) with European- and Chinese reservoir-associated SARS-related CoV strains which 
highlights the sequence diversity in the PLP domain. The alignment shows strains (represented as grey lines) with unique mutations (black 
vertical bars). The ubiquitin-binding surfaces are highlighted in orange (Ub1) and. magenta (Ub2), respectively. The proposed localization of 
Ubi and Ub? is depicted in the cartoon below the alignment. SA- and SR-PLP, which were investigated in the presented study, are highlighted 
in red. Homology modelling for SR-PLP was done to compare the interaction surfaces between SA- and SR-PLP, respectively, and the two 
ubiquitin molecules of a K48-linked di-ubiquitin. Molecular graphics and analyzes were performed with the UCSF Chimera package [71] with 
the crystal structure of SARS-CoV PLP in complex with a K48-linked di-ubiquitin (PDB entry 5e6j) as template and the amino acid sequence of 
SR-PLP as a target. Putative amino acid contact points between SA- and SR-PLP and Ub1 (b) and Ub2 (c) highlight the conformational 
variability of both PLPs in ubiquitin-binding. The tertiary structure of SA-PLP is shown in dark blue, while the SR-PLP tertiary structure is 
shown in lightblue. 


https://doi.org/10.1371/journal.ppat.1007296.g008 





(permit no. CHRPE49/09; A04957) as described previously in [45]. Primary human airway 
epithelial cells were procured from patients who underwent surgical lung resection for any 
pulmonary disease and who gave informed consent. This was done in accordance with the eth- 
ical approval EKSG 11/044, EKSG 11/103 and KEK-BE 302/2015. 
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Cells and viruses 


HEK-293T (Friedemann Weber, University of Gießen) and MA-104 (Friedemann Weber, 
University of Gießen), Vero (Jindrich Cinatl, University of Frankfurt), Vero E6 (ATCC, 
ATCC CRL-1586), Calu-3 (ATCC HTB-55) and RhiLu-hACE2 (provided by author Marcel A. 
Müller [56]) cells were cultivated in Dulbecco’s modified Eagles medium, (DMEM) supple- 
mented with 10% fetal bovine serum (ThermoFisher Scientific), 1% penicillin/streptomycin, 
1% non-essential amino acids, 1% L-glutamine and 1% sodium pyruvate in a 5% CO, atmo- 
sphere at 37°C. HAE were generated and maintained as previously described [64]. RVFV Cl 

13 was a kind gift from Friedemann Weber (University of Gießen). Infection experiments with 
rSCV were done under biosafety level 3 conditions with enhanced respiratory personal protec- 
tion equipment. 


Plasmids 


To ensure high level protein production of SA-, SR- and SO-PLPs in eukaryotic cells the 
codon-usage was optimized based on the human codon-usage frequency. In addition, poten- 
tial splice sites and polyadenylation signal sequences were eliminated before the sequences 
were cloned into the eukaryotic expression plasmid pCAGGS along with a carboxy-terminal 
FLAG epitope tag. Regions of SA-PLP were PCR amplified from parental plasmid 
pPLpro1541-1855 [48]. SR- and SO-PLPs were synthesized by Geneart. Primers, used for clon- 
ing are listed in S1 Table. Site directed mutagenesis was done to change the catalytic cysteines 
into alanines (QuikChange Mutagenesis) and replace the ubiquitin-binding M209 residue to 
arginine (Gibson assembly, NEB) using the indicated primers below (S1 Table). The insertion 
of correct mutations was verified by DNA sequencing. 

Nsp2/3-GFP, pCAGGS-HA-Ub, pISG15-myc and pRL-TK plasmids were kindly provided 
by Ralph S“ Baric (University of North Carolina, USA), Adriano Marchese (Medical College of 
Wisconsin, USA, Min-Jung Kim (Pohang University, Republic of Korea) and Karl-Klaus Con- 
zelmann (University of Munich, Germany), PeDNA3-UbelL, peDNA3-UbcH8 and pcDNA3- 
Herc5 were kind gifts from Robert M“ Krug (University of Texas, USA). P125-Luc, pEF-BOS- 
MDA5-3xFLAG His10, GFP-IRF-3 and pCAGGS were previously described [65]. 


(R)-N-(3-Acetamidobenzyl)-1-(1-(naphthalen-1-yl)ethyl)-piperidine- 
4-carboxamide (3e) 
The synthesis and characterization of SA-PLP inhibitor 3e are described in [51]. 


Generation of SO-PLP sequence 


The SO-PLP sequence was obtained from a Hipposideros bat fecal sample (BUO2-B-F114) as 
previously described by Pfefferle etal. [45]. According to a nucleotide sequence alignment, 
containing related PLP gene sequences, primers were located to the most conserved regions 
within- and downstream of the PLP domain. Two different primer sets were applied in two 
successively performed PCR reactions. The nucleotide sequence information gained was used 
for the design of primers specifically targeting the SO-PLP sequence. The sequencing strategy 
and the primers are shown in S2 and S3 Tables. 


Chimeric SR-PLP-rSCV 


SR-PLP-rSCV was constructed using the previously established cDNA clone [66]. This 
approach was based on bacterial artificial chromosomes (BAC) for keeping the full-length 
CoV cDNA stable. For the construction of the full-length infectious cDNA clone seven 
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subclones (referred to as pAl, pA2, pB, pC, pD, pE and pF) covering the whole SARS-CoV 
genome were generated and assembled in a stepwise procedure. For the generation of SR-PLP- 
rSCV a chimeric SR-PLP subclone, named pBG-ABCD2, containing approximately one-half 
of the SARS-CoV genome and the desired SR-PLP replacement, was used. PBG-ABCD2 was 
joined with subclone pDEF into a full-length BAC cDNA clone and rescued as in [66]. 
SR-PLP-rSCV was sequenced to confirm the presence of SR-PLP and the absence of any fur- 
ther mutations with the following primers: primer for reverse transcription (Brev:5’-TGAA 
CCGCCACGCTGGCTAAACC-3’), sequencing primers (B4622F: 5-CTTAAAGCTCCTG 
CCGTAGTG-3’, BG4792F: 5-TATTAAGGTGTTCACAACTGTAG-3’, BG5631F: 5-AAAT 
TGATGGTGCTCTCTTGAC-3’) 


Virus infection 


Cells were seeded at a concentration of 3.5x10° cells/ml. After 24 h, virus stocks were diluted 
in serum-free medium according to the desired MOI. For virus adsorption 200 ul (24-well) or 
1 ml (6-well) of virus master mix was added to the cells and incubated for 1 h at 37°C. After 1 
h, the virus dilutions were removed and the wells were washed twice with 1x PBS and refilled 
with supplemented DMEM. Supernatants were taken at the indicated,time points and studied 
further. 

For infection of HAE air liquid interface cultures, the apical surface was washed twice with 
200 ul Hank’s balanced salt solution (HBSS) to remove mucus. Virus stocks were diluted in 
HBSS and HAE were infected with an absolute infectious dose of 40,000 PFU. Cells were incu- 
bated for 1.5 h at 37°C in a 5% CO, atmosphere with 95% humidity. After adsorption, virus 
dilutions were removed and the wells were washed three times with 200 ul HBSS. Samples 
were taken at the indicated time points by applying 200 ul HBSS to the apical surface 10 min- 
utes prior tothe actual time points. Basolateral medium was exchanged at 48 hpi. 


Trans-cleavage assay 


Assessment of protease activity was conducted by a trans-cleavage assay. 2x10” cells/ml were 
seeded in a 12-well plate 24 h prior to transfection. 300. ng of PLP-encoding plasmids were 
coexpressed with 25 ng of SARS-CoV nsp2/nsp3-GFP substrate [49]. After 16 h, cells were 
lysed with 100 ul lysis buffer (20 mM Tris [pH 7.5], 150 mM NaCl, 1 mM EGTA, 1 mM 
EDTA, 1% Triton X-100, 12.5mM Na pyrophosphate, 1 mM B-glycerophosphate, 1 mM Na 
ortho-vanadate, 1 mg/ml leupeptin, 1 mM PMSF) and lysates were separated by SDS-PAGE 
using a semi-dry transfer blotter. After protein transfer, the membrane was blocked by 5% 
dried milk in TBST buffer (0.9% NaCl, 10 mM Tris-HCl [pH 7.5], 0.1% Tween 20) for 1 h at 
room temperature. The membrane was incubated with a rabbit-antibody directed against the 
GFP-epitope tag (ThermoFisher Scientific). The membrane was washed three times for 15 min 
in TBST buffer. Next, the membrane was incubated with HRP-coupled donkey anti-rabbit sec- 
ondary antibody (SouthernBiotech). After 1 h, the membrane was washed three times for 15 
min in TBST buffer. To confirm PLP expression, the membrane was probed with mouse anti- 
FLAG (Sigma-Aldrich) followed by goat anti-mouse HRP-coupled (SouthernBiotech) antibod- 
ies. Mouse anti-B-actin (Sigma-Aldrich) was used to detect host cell proteins as a loading control. 
Secondary detection was performed using goat anti-mouse HRP antibody (SouthernBiotech). 


Biosensor assay 


The assay was performed using HEK-293T cells in black 96-well plates with clear bottom. Cells 
were transfected with 37.5 ng pGlo-30F-RLKGG [50] and 50 ng PLP-expressing plasmids 
using Fugene HD (Promega) according to the manufacturer’s instructions. At 14 hours post 
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transfection (hpt), cells were equilibrated with GloSensor reagent (Promega). For inhibitor 
quantification, cells were equilibrated at 13 hpt with GloSensor reagent as indicated above and 
at 14 hpt diluted inhibitor or DMSO was added. Luminescence was measured after 1 h for the 
following 6 h. For data analysis, the fold luciferase induction was calculated independently for 
each PLP by calibrating the values, obtained for every time point to the respective starting 
values. 


DUB activity assay 


HEK-293T cells (2x10 cells/ml) were transfected in the 12-well format with 300 ng of HA-ubi- 
quitin (HA-Ub) and 50, 100 or 200 ng PLP-expressing plasmids using TransIT-LT1 transfec- 
tion reagent (Mirus). At 18 hpt, cells were treated with 100 ul of lysis buffer. Proteins were 
separated by SDS-PAGE and analyzed by Western blotting as described above. Western blot 
analysis was done using mouse anti-HA serum (Covance) and goat anti-mouse HRP (South- 
ernBiotech) antibodies. 


DeISGylation assay 


HEK-293T cells were transfected with plasmids encoding for myc-epitope tagged ISG15 (250 
ng of pISG15-myc) and its conjugation machinery, comprising a set of three ligases (UbcH8/ 
125 ng pUbcH8, Ube1L/125 ng pUbelL and Herce5/125 ng pHerc5) as in [40]. PLP-encoding 
plasmids were cotransfected in amounts of 100 ng. After 18 h, cells were lysed and separated 
by SDS-PAGE as indicated above. The extent of ISGylated cellular proteins was analyzed by 
Western blotting using mouse anti-myc (MBL Life science) and goat anti-mouse HRP (South- 
ernBiotech) antibodies. 


Dual luciferase assay 


2x10° HEK-293T cells/ml were transfected with plasmids encoding Renilla (RL) and Firefly 
(FF) luciferases, The FF luciferase gene was under control of an IFN-B promoter (p125-luc). 
RL luciferase gene was cloned behind a herpes simplex thymidine-kinase (TK) promoter. IFN- 
B promoter activation was triggered by the overexpression of the cellular MDAS5, which led to 
an auto-activation of the IFN pathway [67]. Transfection of DNA plasmids was done with 
Fugene HD (Promega) according to the manufacturer’s instructions. Plasmid amounts are 
given for the 24-well plate format and are as follows: pRL-TK (5 ng), p125-luc (250 ng), pEF- 
BOS-MDA5-3xFLAG His10 (100 ng).and 1 to 50 ng of PLP-expressing plasmids. Cells were 
lysed at 17 hpt. 20 uLof luciferase-containing lysate was transferred to an opaque, white 
96-well microtiter plate and used for assessment of luciferases activity by a bioluminescence 
detection reader. 


Translocation assay 


The IRF-3 translocation assay was performed as described elsewhere [65]. Briefly, MA-104 
cells were transfected with plasmids encoding 250 ng GFP-IRF-3 and either EV or 1000 ng of 
FLAG-tagged PLP-expressing plasmids using Fugene HD (Promega). To activate the IFN 
pathway, cells were infected with RVFV Cl 13 and an MOT of 5 at 17 hpt. At 8 hpi cells were 
fixed with 4% paraformaldehyde and permeabilized with 0.1% TritonX-100. Immunofluores- 
cence analysis was done as in [68]. Samples transfected with EV were treated with anti-RVFV 
mouse serum (Friedemann Weber, University of Gießen, [53]) and goat anti-mouse Cy3 (Dia- 
nova GmbH) secondary antibody. PLP-expressing cells were stained with mouse anti-FLAG 
(Sigma-Aldrich) and goat anti-mouse Cy3 (Dianova GmbH) antibodies. Samples were 
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analyzed by a fluorescence microscope (Zeiss). Depending on the number of transfected cells, 
at least three images were taken, and the number of cells double-positive for GFP and FLAG 
was divided by the number of cells showing IRF-3 nuclear translocations. 


Plaque titration assay 


3.5x10° VeroE6 cells/ml were seeded in a 24-well plate 24 h prior to infection. A 1:10 serial 
dilution of samples was generated, 200 ul dilution was added to the cells and incubated for 1 h 
at 37°C for adsorption. After virus samples were removed, 500 pthoverlay (2.4% avicel diluted 
in 2x DMEM) was added to each well. The overlay was discarded at 3.days post infection (dpi) 
and cells were fixed for 30 min in 6% formaldehyde. The cells were washed once with 1x PBS 
and stained with crystal violet working solution for 15 min. Plaque forming units were deter- 
mined from at least two dilutions for which distinct plaques were detectable. 


Real-time RT-PCR assays 


Viral RNA was extracted with the NucleoSpin RNA virus isolation kit (Macherey-Nagel) after 
the manufacturer’s instructions. 1 ul of viral RNA was applied to each reaction. Quantification 
of genomic SARS-CoV RNA was done using the SuperScript III one-step reverse transcrip- 
tase-PCR system (Invitrogen) with the Platinum Taq DNA polymerase according to the manu- 
facturer’s recommendations and these primers (SARS-F: 5’-GCCGCGAAGAAGCTATTCG- 
3’, SARS-P: Fam-5’-ACGTTCGTGCGTGGATTGGCTTTG-3’-BHQ, SARS-R: 5-AGTTGC 
ATGACAGCCCTCTACA-3’). An.established SARS-CoV standard curve was used in each 
run [69] to quantify RNA copies per ml. 

To analyze IFN-B expression in Calu-3 and Rhinolophus alcyone cells (provided by authors 
Christian Drosten and Marcel A. Müller), mRNA was extracted with the NucleoSpin RNA iso- 
lation kit (Macherey-Nagel). Real-time PCR was done as described previously [70] using the 
following primers: Hu-IFNB1-F: 5’-AGGATTCTGCATTACCTGAAGG-3’, Hu-IFNB1-P: 
Fam-5’- TCCACTCTGACTATGGTCCAGGGA-3’-ZEN and Hu-IFNB1-R: 5-GGCTAGGA 
GATCTTCAGTTTCG-3’ (Calu-3) and Rhi-IFNB1-F: 5’°-AAATAATGGAGGAGGAAAACT 
TCAC-3’, Rhi-IFNB1-P: Fam-5’-CGAAACATGAGCACGCTGCACCTG-3’-BHQ and Rhi- 
IFNB1-R:5’- C@CCTGATCCTTAGGTAGTAATTCT-3’ (Rhinolophus alcyone). To deter- 
mine the induction of IFN-f relative to TBP (Calu-3) and f-actin (Rhinolophus alcyone) the 
2-4“ method was applied [55]: 


Alignments, phylogenetic analysis, similarity/identity calculations 


Phylogenetic analysis was performed using the neighbor joining algorithm in Geneious 9.1.5 
under the assumption of a Tamura-Nei genetic distance model. Alignments were done based 
on the amino acid codes by the BLOSUM62 algorithm in the Geneious 9.1.5 software package. 
The identities and similarities of the listed proteins were calculated from the amino acid align- 
ment using the BLOSUM62 substitution matrix with a threshold of 1. 


Accession numbers 
SO-PLP: MG916963 


Supporting information 


S1 Fig. Analysis of SR-PLP splicing products by Western Blot analysis. HEK-293T cells 
were transfected with a plasmid expressing non-optimized SR-PLP. For detection of potential 
splicing products the amino-terminus contained a Flag- and the carboxy-terminus contained 
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an HA-tag. Protein expression was confirmed by Western Blot analysis using Flag- or HA-spe- 
cific antibodies produced in mouse, respectively. Secondary detection was done using anti- 
mouse antibodies coupled with horseradish peroxidase, which were produced in goat. 

(TIF) 


S2 Fig. Comparative inhibition of rSCV and chimeric SR-PLP-rSCV by compound 3e. The 
number of viable cells after inhibitor treatment was determined by a luminescence-based cell 
viability assay. Therefore, Vero cells were seeded in an opaque 96-well plate and after 20 h, 
cells were incubated with either DMSO or a serially diluted compound 3e using DMEM as dil- 
uent. After 24 h, 100 ul of detection reagent was added and incubated.for another 15 min. 
Emitted luminescence was recorded by a detection reader with an integration time of 1s. 

(TIF) 


S1 Table. Primers for cloning of the constructs used for heterologous expression analyzes. 
*Underlined nucleotides were added for cloning purposes. "FLAG tag is highlighted in bold. 
For: Forward primer; Rev: Reverse primer. 

(DOCX) 


S2 Table. Amplification strategy for the generation of SO-PLP sequence from bat fecal 
sample BUO2-B-F114. SSIII one step: SuperScript II] one-step RT-PCR system with Platinum 
Taq DNA polymerase. 

(DOCX) 


S3 Table. Primers for the amplification and sequencing of SO-PLP. For: Forward primer; 
Rev: Reverse primer. The single-letter code system for degenerated bases established by the 
international union of pure and applied chemistry (IUPAC) was used for definition [72]. 
(DOCX) 
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